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1. INTRODUCTION 

Nowadays, efforts by governments and businesses are further promoting to propel the mass adoption 
of renewable energy in order to steer the production of electricity from clean energy sources. Among these 
energies, the renewable wind energies are becoming used widely in this context, because it is cost-effective, 
sustainable, clean, and does not produce greenhouse gases. 

Recently, many wind turbine technologies have been developed, where the main goal is to ensure a 
good quality of electrical power, and an optimization of the energy resource that is extracted from the wind. 
Therefore, the choice of the most appropriate electrical generator for a wind turbine is an important task. 
Appropriately, the most commonly used generators in wind power generation, are based on a permanent 
magnet synchronous generator (PMSG), for its many advantages compare to other machines [1]. The PMSG 
does not require electrical excitation, as the magnetic field is produced by permanent magnets, rather than by 
the coil. Hence, PMSG’s does not need slip rings nor brushes, which reduces the weight, implementation 
costs, maintenance and no field copper losses. Furthermore, there is the possibility to avoid gearbox 
connection to the turbine. It will result in high dynamic performance, and high power with a wide operating 
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speed range [2], [3]. As is well known, the complex dynamic behavior describing wind turbine systems, 
suggests robust control strategies to achieve the desired performance and ensure high stability [4], [5]. 

In previous findings, many control methods are proposed including linear controllers that are limited 
in the ability, to achieve acceptable control performances. Nonlinear controllers are better adapted to the 
wind energy conversion system (PMSG-WECS) under different wind speed variations [6]. Makhad et al. [7] 
design a novel integral backstepping strategy, to optimize wind energy in the event of abrupt changes in wind 
speed. Dahbi et al. [8] clarified the wind turbine system by combining maximum power point tracking 
(MPPT)-pitch angle control, using artificial neural networks to enable the network to be supplied with stable 
and efficient electrical energy. In the power system, the wind farm is more sensitive to disturbances. In this 
context, Chen et al. [9] designed MPPT control scheme based on nonlinear adaptive control to estimate the 
lumped perturbation term according to a high-gain perturbation observer. Moradi et al. [10] design a Hx 
controller to improve the wind turbine power with robust stability and a less oscillatory behavior compared 
to the classical PID controller. 

Moreover, fuzzy logic is well suited to modelling a nonlinear system according to [11]-[14]. In this 
context, Allouche et al. [15] designed a Takagi- Sugeno reference fuzzy model, in order to generate the 
optimal trajectory corresponding to only the maximum power. However, in that work, while also operating 
under time varying wind speed, the case of strong wind is not considered, despite the fact that the T-S fuzzy 
control can provide a wide large of control gain variation. 

Based on the principles of the T-S fuzzy controller, which can provide an effective representation of 
complex and nonlinear systems, involving system uncertainties and external disturbances, the suggested 
solution in this work is to attenuate influence of external disturbances. In this article the tracking control 
design based on the augmented T-S fuzzy model by defining a lumped perturbation, is described using Hə 
performance where the tracking error must be equal or less than a prescribed attenuation level. In this 
approach, the direct-drive permanent magnet wind turbine system is described by an aggregation of linear 
models, interconnected through membership functions. Then, local linear compensators for each subsystem 
to achieve the desired objective determine the global controller system. 

Under a wide range of the wind speed, two different operations must be considered. The first one is 
a maximum power point tracking (MPPT) controller, established to optimize the generator speed. In this way, 
the captured energy is set at its maximum level, under rated wind speed. On the other hand, pitch angle 
operation control strategy is focused in limiting the turbine output power and the generator speed 
respectively, exceeding the wind speed rated value. Furthermore, in our approach, the model based T-S fuzzy 
controller is designed via parallel-distributed compensation (PDC), to determine the state feedback fuzzy 
controller, from the proposed augmented T-S fuzzy models. So, to improve the suggested controller, as a 
contribution, a stability criterion is derived from Lyapunov’s direct method and, the H» tracking control 
performance can then be guaranteed at the same time. In this context, the feedback gains are solved very 
efficiently by convex optimization linear matrix inequalities (LMIs). Also, in order to provide optimal 
dynamic performance in term of convergence and robustness, the proposed method takes also the system 
parameter uncertainties of the mathematical model into account during the design process as another 
significant contribution. 

In the following, section 2 includes a T-S fuzzy logic description for the wind turbine model based 
on the PMSG. In section 3, a tracking augmented T-S fuzzy controller is designed via the parallel distributed 
compensation (PDC) scheme and the Lyapunov stability criterion is presented for the closed-loop system to 
provide efficient power conversion. Then the robust tracking control performance Hə is described, to 
guarantee the stability of the system and attenuate the disturbances, respectively. Linear matrix inequality 
(LMI) technique is adopted to solve for the feedback gains of the fuzzy controller. Section 4, presents 
simulation results applied for a 5 MW permanent magnet synchronous generator. The conclusion and future 
scope close the article. 


2. WIND ENERGY CONVERSION SYSTEM 
2.1. Wind energy model 
The mechanical power harnessed by a horizontal-axis wind turbine can be expressed in (1) [15]: 


1 
Pn = = Cp (A, B)pmR ey (1) 
where R, is the rotor radius, p is the air density, Vw is the wind speed and C,(A,B) is a power coefficient 


presented as a nonlinear function that depends on tip speed ratio A and the pitch angle of the blades B as 
shown in Figure 1. 
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The empirical and nonlinear equation of C(A, fp) based on the turbine characteristics are given by 
the expression in (2) [16]: 
— 151 2.14 -124 
Cp(A, B) = 0.73 (Œ — 0.586 — 0.026?"*p — 13.2) e~a (2) 


1 
A = 1 0.003 
A-0.02B8 B341 


The tip speed ratio is defined as (3): 


— NQmRt 
A= Sn (3) 
where (Qm 1s the turbine speed. 
The typical variations of C(A, P) from (2) with the tip-speed ratio À for various values of pitch 
angle 6 are illustrated in Figure 1. 
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Figure 1. Coefficient of performance as a function of Tip Speed Ratio 


The PMSG speed control based on MPPT (maximum power point tracking) control consists in 
maximizing the power converted. Thus, the optimal speed of the turbine must be adjusted as (4) [17], [18]: 





hoptYw 
Qmopt = Ri (4) 


for each wind speed with maximum C, = 0.44 and optimal à = 6.89 as shown in Figure 1. 
The purpose of the MPPT can therefore be provided by setting the electromagnetic power extracted 
from the wind at the maximum value given by (5) [15]: 


5 

_ 1pTRECymax n3 

Prax — 2 43 Qa (5) 
opt 


When the wind velocities are higher than rated, the blades are rotated by a control device based on pitch 
control as shown in Figure 2, to reduce the aerodynamic power captured by the wind and maintain the output 
power of PMSG at its rated value [19], [20]. 
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Figure 2. Pitch angle controller 


2.2. PMSG model 
For a non-salient PM machine, the dynamic model of the PMSG with rotational topology is given in 
the d—q synchronous frame using the Park’s transformation as defined in (6) [16]: 


Va = Rig + La <4 — Welgi, = 


: di : 
Vq = Riq + La = + We(Laia + $y) 


Where Vy, Vj are the stator voltages in the d-q axis, lg, iq are the currents in the d-q axis, We = ply is the 
electrical rotation speed, R is the stator resistance, ®, the flux linkage of permanent magnets and p is the 
number of pole pairs. The dynamic equation of the wind turbine is given by (7) [16]: 


dQm 


dt Sg lag Jin (7) 


J is equal to the summation of inertia moments of the turbine and the generator, f is the friction coefficient, 
T, and Tm prensents respetively the electromagnetic and aerodynamic torque. 
The equation of the electromagnetic torque produced by the machine is defined as (8): 


3 r 
Te = 5 PP ela (8) 


Using (6) and (7), the dynamic model of the PMSG-WT in d-q reference frame can be described in 

the following nonlinear state space form (9) [24]: 
x(t) = f (x(t), ul) + g(x®)w@ 

(9) 

y(t) = q(x) 


where x E R” is the state vector, u E R™ is the control input vector, w(t) denotes the external disturbance, 
y E R” is the mesured output, f(x), g(x) and q(x)are nonlinear with appropriate dimensions. 

In order to express the nonlinear model of the machine as a T-S fuzzy model, with the measurable parameters 
as decision variables, (6) is described in the following nonlinear state space form (10) [24], [25]: 


f (t) = A(O,,) x(t) + B u(t) + E w(t) (10) 
y(t) =C x(t) 
where 
x(t) = [ig(t) ig(t) mE) ed 
= DQM 0 = 0 0 
A(Qm) = —pQ.n = ae ,B = 7 ij, E= a ,C=[0 0 1],ut)= o 
i f 0 a J ' 
2 J J 


(12) 
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2.3. Takagi-Sugeno (T-S) fuzzy model 

The T-S fuzzy model of PMSG is established based on the dynamic model given by (9). Thus, for 
ensuring the control of the stator converter side generator, the system is approximated by the T-S fuzzy 
dynamic model composed of r rules in (11). 

The it? rule of the fuzzy model is as (13) [22], [23]: 


IF v, (t)is M,; and ...and v; (t) is M THEN x(t) = A; x(t) + B; u(t) + Ew(t) (13) 


i= 1,2,...,r; Mj; denotes the fuzzy sets and represent the grade of membership of v4 (t) in Mj;. 
r is the number of IF .... THEN rule and v;(t) = (v4, v2, ..., vj) E R/are the premise variables [24] 
A; B; are the local subsystem matrices. 

Using the singleton fuzzifier, the resulting overall fuzzy system is inferred as (14) and (15) 


x(t) = Mini UUE) Axl) + Bru) + E,w(t)) 
r (14) 
y(t) = din (UC) Cx) 
l — _ Mii) 
HO) = FO = 
Furthermore, it should be noted that the expression u;(v(t)) for all i satisfy (16): 
dizi Hi(u(t)) = 0 (16) 


With Xf 4;(v(t)) = 1 for all t > 0 where y,(v(t)) = 0 for i = 1,2, ...,r 
Considering the rotor speed chosen as premise variable, two rules are deduced and the sub-matrices can be 
written as (17): 


R R 


=, L, PQmax 0 = La PQmin 0 
d d 
is aR i 6 a aR cs BEE 
A, = PQ max Ey La , A> = PQ min Le Ly ) (17) 
0 3 PPF = f 0 SPs = f 
2 J J 2 J J 
=. 6 


La 
Bı = Bz =19 1|, E = E, = 


Lq 
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| 
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Figure 3 illustrates clearly the principle of the proposed fuzzy controller which corresponds to the 
maximum energy captured from the wind and its limitation at the nominal value. 
In this regard, to simplify the tracking controller design, a set of new optimal trajectory is (18): 


Xot) = EEA loot C) Umont (t)]" (18) 

x(t) — Xopt (t) > Oast > 0 

X(t) T x(t) = Xopt (t) 

X(t) is defined as the tracking error and its time derivative is given by (19): 

X(t) = X(t) — žopt (t) = Dh iVE) (AZC) + Bult) + Eiwt))—žopt(t) (19) 
The new fuzzy controller t(t)is substituted in (14) is defined by(20): 


ray ilU) BTE) = Er- ui lult) )(AiXopt(t) + Byw(t)) — žopt (t) (20) 


(19) is rewritten in the following compact form (21): 


X(t) = Via mi (v(t) ) (AZE) + Bit(t) + E;w()) (21) 
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Figure 3. Generator-side control scheme 


3. ROBUST AUGMENTED T-S FUZZY CONTROLLER DESIGN 
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In the next section, this controller is designed following the Parallel distributed compensation (PDC) 


technique employed for the T-S fuzzy model (21) [25], [27]. 
Hence, the it” rule fuzzy controller can be described as (22): 


Rt: IF v,(t)is M4; and ...and v; (t) is Mj,THEN t(t) = —K;X(t), fori = 1,2,. 


aro (22) 


Where K; € IR"*™ is the local state feedback vector. Then, the overall fuzzy PDC controller (22) is inferred 


as (23) [26]: 
T(t) = — Ye Mi(v(t)) KZE) 
Substituting (23) into (21) yields the closed-loop fuzzy system as (24): 
X(t) = Vie Dj-1 Hi(vO©) wj(v@) (Ai — BKZ) + Eiw() 
let us denote the integral state error vector as (25): 
O) = f x(t)dt 
The augmented system containing the new fuzzy controller is presented as (26): 
X(t) = Dir mi(v()) (AX O + BTO + E,W) 
Where the new feedback controller is inferred as (27): 
T (t) = — Vie. mi (v()) Ki X(t) — LE, wi (v(t) Fe (0) 
Thus, (27) is rewritten as (28): 


X(t) 


T(t) = — Xi- mi (v(o))[K; Fi] P 
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(26) 


(27) 
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7 J4 0] g -Eq p -= [El m = T 
A =F o= lah E= [il 70 = wo o (28) 
Substituting (28) into (26) yields the augmented closed-loop continuous T-S fuzzy model as (29): 


Š O = Dh Vier MVE) WA) (GX O + WO) (29) 


3.1. H» tracking control performance 

The optimality criterion by Hæ control performance consists to synthesize a controller such that the 
equilibrium of the closed-loop system (29) is stable and the attenuation of external disturbances is guaranteed 
considering the Hœ performance as (30) [25], [27]: 


Jt XT dt X < X70) X (0) + 8? LW" (Wejdt (30) 
Where t; denotes the final times and ô is a specified disturbance constant. 
In order to ensure the asymptotic stability of the augmented closed-loop model (29) guaranteeing the 
Ho tracking performance described (30) for all disturbance, the gains K,is obtained from Gy. 
Lemma. 

The augmented closed-loop system described by (29) is asymptotically stable with the prescribed 
Ho performance inequality (30), if and only if there exists a matrix P = PT > 0 


such that (31) 


A, TX +X A," — B,M; — M,"B," E; X 


p _s27 9 | <0 (31) 


X 0 =] 


Proof. 
Let the Lyapunov function for the augmented closed-loop system be defined as (32) 


V =XT(t)PX(t)>0 (32) 
Taking the derivative V(X (t)) will be required to satisfy (33) [23], [26] 
XT (t) PX (t) + XT PRT (t) <0 (33) 
H output tracking index 6 is achieved for the resulting augmented closed-loop fuzzy if system (23) holds 
P(O) <- OEO ++ CWO" WO <0 (34) 
Once the following matrix inequalities (for all 1) are satisfied, 1.e. 


Whar Mi(v(t)) (Gyf (O + EW E)) PX O + X (OP (GX O + EWE) +T (OX - 
6*W(t)? W(t) < 0 (35) 


r ulv) (OLG,P + P GJIT (t) + W(e)[E,PIX (t) + P EWA) — 5?2W(t)? W(t) <0 








(36) 
FOPO ae GP wes +1 P Yi ui (u(t))E; XO] eo = 
"1 ui (u(t) )E; — §7] W(t) 
ae (G.P+PG,)+1 P Yt. u,(v(t))E, Ai 
Dros Mi(v(t))E; P =O 





Augmented robust T-S fuzzy control based PMSG wind turbine improved ... (Naoual Tidjani) 


592 o ISSN: 2088-8694 


Hence, it suffices to check that 





r r 
—— —T — 
> uilult)) (GP +PG ) P > wi; (u(t) )E; 
a iat I 0 (39) 
: E + p | <0 
> ui(v(t))E; P — 6°] 
j=i 
r r 
SE —T — 
ulu) CGP +PG ) P > ui(ve))E; (40) 
ie a =e +[]U 0] <0 
2 mi(v(e)) EP -ôI 
j=1 
T r 
— —T = 
u(v(t)) (GP+PG,) P > ui(ve))E; I 
j=1 j=1 (41) 
r , <0 
» ui (u(t) )E; P —§?] 0 
j=l 
I 0 =] 
GP+PG PE, 1 Á 
EP -8I o0|<° 
I 0 -I 
-1 9 ol [GP+PG" PE 1][p-2 0 0 a 
o I o EP = —6?] o || 0 1 oj<0 
I 00 I 0 —] ime 
Thus, (43) 1s equivalent to (44) 
P-1(GP+PG,)P-! P-1PE, P~! 
BE," PP- 671 0 i 
pot 0 z. 


Let's consider X = P~1 and M; = K,P™+ ; the control gains are obtaining by solving the LMI obtained by (45) 


XA +ĀX-B MoM Be Fo X 
E -6I 0 
X 0 el 


(45) 


3.2. Virtual disired variables system 
By using the fact that g(t) = X}=1 ulv (t)) (20) is rewritten in the compact form (46) [25]: 


g(t)Byr(t) = g(t) (AiXope(t) + BruCt)) — xope(t) B(x(t)) (a(t) — u(t) = A(x(t))xope(t) + 
Xope(t) (46) 


where 


Ai(x(t)) = g(t)A;, Bi(x(t)) = g@B; (47) 
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R 


1 —— po 0 
La j “a i D igopt (t) ldopt (t) 
DYE ; . 
0 ~ (u(t) E t(t)) a Pem 7 Iy E Pa igopt (t) T tqopt (t) (48) 
0 0 0 See =F Omopt (t) Omopt (t) 
2 J J 


In order to provide the maximum torque per ampere i dopt 1$ always kept at zero. And from the measured 
optimal speed, Xopt(t) can be deduced from (48). 
Therefore, the desired state variables for the tracking control are set to be 
0 
2] (¢ f 
Xopt t) = pof (Can (t) + 7 Smopt (t)) (49) 
nopi (t) 


The required components of the rectified voltage vector are derived from the integral controllers for each d-q 
current component. Consequently, from (6) the machine voltages are described as (50): 


| Va = —PpLamiqopt + Rigopt + idont + Tq (50) 
Vq = PP mope + Riqopt + LaPmtaope + Tq 


4. RESULTS AND DISCUSSION 
Matlab simpower system environment is employed to validate the elaborated model and control of 
the conversion system wind energy. The PMSG and turbine parameters are listed in Table 1. 


Table 1. Specification of PMSG and turbine [14] 


PMSG parameters Values Wind turbine parameters Values 
Generator rated power 5 MW Turbine blade radius 58m 
Generator rated voltage 3.3 kV Optimal tip speed 6.89 
Generator rated current 1.515 kA Maximum C, value 0.44 
Pole par number 75 The air density 1.225 kg/m?* 
Generator stator resistance 6.23 e-003 Q Rated wind speed 12.12 m/s 
Generator d — q axis inductance 4.229 e-003 H 
Permanent magnet flux 11.1464 Wb 
Generator inertia 2 e+5 Kg. m? 


The evolution of wind speed is illustrated in Figure 4, which shows that the speed varies around 
10m/s + 30% during the time interval 0 to 60 s. 

Depending on the wind velocity variation, two operations can be triggered according to the 
appropriate control strategy. 
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Figure 4. Wind Profile Figure 5. Generator rotor speed 
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Figure 5 shows that below the rated wind speed, the generator speed is adjusted to be optimal when 
regulating the system under the specification of maximum power, and is limited at the rated value beyond the 
rated wind speed without overshooting. This shows the robustness of the method against the wind speed 
disturbance. 

Figures 6(a) and 6(b) show that below the rated wind speed, the energy captured from the wind is 
maximized from MPPT operation, by adjusting the optimal rotor speed that corresponds to the maximum 
power. In this context, the power coefficient is kept constant at its maximum value at about 0.44, and the 
pitch angle B is preserved at zero for each MPPT operation, despite the variation of the wind. 

Whereas beyond the rated wind speed value V, = 12.12 m/s, the turbine is protected against 
damage and the blades are rotated by adjusting the pitch angle. As a result, B increases with subsequent 
degradation of tip speed ratio A in order to reduce power coefficient Cp to a value less than 0.44, as shown at 
pitch control operation. 
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Figure 6. (a) Power coefficient, (b) Tip speed ratio and pitch angle variation 


Figure 7 shows current responses in d q frame, as the currents track well the reference trajectory. In 
this context, the quadrature current i, (Figure 7 (a)) varies according to the generator speed, taking into 
account the two wind operating modes (MPPT and pitch angle). 

The direct current id (Figure 7 (b)) tracks well the optimal reference (at zero) in order to provide the 
maximum torque per ampere. 
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Figure 7. Stator current in d-q frame, (a) Quadrature stator current, (b) Direct stator current 


Figure 8 depicts the stator currents of the PMSG in sinusoidal form. It can be noticed that stator 
currents have a similar shape as the wind speed during the time intervals [Os, 5s] and also is limited, beyond 
the nominal wind speed as it shown during the interval [5s, 8s]. 
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In order to guarantee the control performance of the closed-loop T-S fuzzy system, the model of the 
system takes into account the variation of the machine parameters (R varies by 30 %, Lg = Lq varies by 
20%), that change at the same time during the time simulation. 

Figure 9 shows that despite variations in the PMSG parameters, once the rated power of the 
generator is reached at the nominal wind speed, the robust T-S fuzzy controller track the output power at its 
nominal value (5 MW) accurately and rapidly even though there is uncertainty. This indicates robustness 
against parameter uncertainties as desired in the problem discussion. The control gains obtained are presented 
in Table 2. 
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Table 2. Feedback gains of the augmented T-S fuzzy controller 


¢ = [1080 0 AN M —0.1574 a 
i D 1080 7480.3 1 199309 53.5681 481.1057 
a =| 90 616. | i =| 44.9273 398.997 


5. CONCLUSION 

This work presented a contribution to the Takagi-Sugeno fuzzy control, designed to control the 
variable-speed wind energy conversion system based on the PMSG. In this context, an MPPT technique was 
applied in order to make the system operate with maximum efficiency. Also, a pitch angle control strategy is 
proposed to avoid damage to the machine in strong winds so as to limit the power of the wind generator at its 
nominal value (SMW). The developed robust augmented T-S fuzzy control was designed via parallel- 
distributed compensation (PDC) involving system uncertainties and external disturbances. Moreover, the 
H tracking control technique performance design is characterized to reject unknown disturbances. The 
stability of T-S fuzzy system is formulated into LMI techniques from a Lyapunov approach. The desired 
virtual variable is introduced to benefit the control design. The very high-performance results of numerical 
simulations, reveals the validity of the proposed method control in terms of robustness and saturation of the 
power. As a perspective, the observer-based finite time H.. control will be integrated. 
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